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Abstraci0; We show that certain integral positive definite symmetric 
tridiagonal matrices of determinant n are in one to one correspondence with 
elements of (Z/nZ)*. We study some properties of this correspondence. 

In a somewhat unrelated second part we discuss a construction which 
associates a sequence of integral polytopes to every integral symmetric ma- 
trix. 



1 Introduction 

Let ai, . . . , a„ be a finite sequence of elements in a commutative ring. We 
denote by T{ai, . . . , On) tlie symmetric tridiagonal matrix 
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of size n X n with diagonal entries ai, . . . , a„ and sub and super diagonals 
consisting only of ones. 

Leighton and Newman discovered in [2] the remarkable fact that the 
number of such matrices of size n x n which are unimodular, integral and 
positive definite is given by the n— th Catalan number ( ^) / (I + n), see also 
Exercice 6.19 nnn in [5j and [4J. One should mention that this result was 
already implicitely contained in [Tj, see the solution to problem (18). 

In the present paper, we are interested in the number of such matrices 
of arbitrary size and determinant which are integral and have diagonal 
entries in {2,3,4, . . . }. Such matrices are always positive definite, the case 
ai = a2 = ■ • • = On = 2 defining the root lattice of type An. 



^Keywords: Tridiagonal matrix, symmetric positive definite matrix, Catalan number, 
Euler totient function, Euclid's algorithm, SI/2(Z), root system, integral polytope, finite 
simple graph. Math, class: 11A05, 11A55, 11AC20, 15A48,52B20 
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Theorem 1.1. The number of finite symmetric tridiagonal integral matrices 
of the form T(ai, 02, • • • ) with diagonal entries > 2 and determinant N is 
given by Euler's totient function 



(where the product is over all prime divisors of N) counting the number 
tt((Z/A'"Z)*) of invertible elements of the ring TLjNTL. 

The proof is bijective: Given (a, A^) for a in {1, 2, . . . , — 1} invertible 
modulo A^, we construct a finite sequence (ai, • • • ) giving rise to a matrix 
r(ai, Q!2) • • • ) with the desired properties. This construction is based on 
properties of SL2(Z) related to Euclid's algorithm. 

The organisation of the paper is as follows: 

The next section exposes a few useful (and probably well-known) facts 
concerning SL2(Z). 

The proof of Theorem 11.11 is contained in Section [3l 

Section H] discusses some algorithmic aspects of the construction wich 
associates to an element a £ (Z/NZ)* the corresponding tridiagonal matrix 
of the form T(ai, . . . ). 

Section [5] describes properties of the function a 1 — > tr(T) — dim(T) for 
a € {Z/NZy. 

All matrices described by Theorem 11.11 are definite positive. Section [6] 
discusses briefiy Euclidean lattices admitting the matrices of Theorem 11.11 
as Gram matrices. 

The final Section [7| is not closely related to the rest of the paper. It 
describes a construction of a sequence of integral polytopes associated to a 
symmetric integral matrix and discusses some easy properties of this con- 
struction. The first polytope obtained by applying this construction to tridi- 
agonal matrices with I's on the sub and super diagonal end zeros elsewhere 
yields for instance integral polytopes with vertices enumerated by Catalan 
numbers. These polytopes are however distinct from associahedrons (or 
Stasheff polytopes). 

2 SL2 and the determinant of T{ai, . . . , a„) 

We consider the map 





extended to finite sequences by setting 



M(ai, 02, ...,«„) = M{ai)M{a2) ■ ■ ■ M(a„) . 
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Denoting by 

|r(ai, . . . ,a„)| = det(T(ai, . . . ,a„)) 
the determinant of T{ai, . . . , a„) we have the fohowing result: 
Proposition 2.1. We have for all n>2 the identity 



M(ai, . . . ,a„ 



-|r(a2, 
|r(ai, . 



,an-i)| 



-|T(a2, 
|r(ai, . 



• ,an)| 



where we use the convention \T{a2, ■ ■ ■ , ctn~i) \ = 1 for n = 2. 

The proof is a straightforward computation left to the reader. 
Proposition 12.11 implies easily the identity 

|r((ai, . . . ,an+i)\ \T{a2, ■ ■ ■ ,an-i)| = \T{ai, . . . ,a„)| \T{a2, ■ ■ ■ ,an+i)\- 1 

which is a particular case of Dodgon's condensation formula. 
Another useful result is the following "symmetry" . 



Proposition 2.2. // 



M{ai,a2, . . . , Or 



then 



M{ar. 



,Oi2,aij 



Proof Introducing the involution V = 
shows VM{a)V = M{a)-^. This yields 

VM{ai,...,an)'W = M{a 
which establishes the result by computing 



d 

1 

1 



a short computation 
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1 
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1 

1 



□ 



The easy identity 

M{x)M{y) = M{x + l)M{l)M{y + 1) 



-y 

xy — 1 



shows that the set of integral positive definite matrices of the form T(ai, a2, ■ ■ ■ 
with given determinant is infinite without further conditions. 

Leighton and Newman consider matrices of the form T(ai, . . . , a„) which 
are integral, positive definite and have fixed size n x n in order to get finite 
sets. In the present paper, the size of the matrices T{ai,a2, . . . ) is arbitrary 
but the diagonal coefficients ai, 02, . . . are restricted to the set {2, 3, 4, ... } 
of natural integers > 2. This restriction is motivated by the following result. 
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Proposition 2.3. Let ai,a2, ■ ■ ■ ,ai be a sequence of I > 1 real numbers 
cii>2 for i = 1, . . . ,1. Consider 



M{ai,a2, ■ ■ ■ 



,ai) 



( 



—a 



c 



d 



b 



) 



We have then the inequalities 



< a < min(6, c), max(6, c) < d and c — a < d — b . 



Corollary 2.4. We have 



det(r(ai,a2, • • • ,az_i)) < det(r(Qi, 02, • • • a;)) 



if ai, . . . ,ai are real numbers > 2. 

In particular, the number of matrices of the form T{ai,a2, ■ ■ ■) with 
determinant N and integral diagonal coefficients ai, 02, • • • C {2, 3, 4, ... } 
is finite. 

Proof The first part follows by applying Proposition 12.11 and the in- 
equality c < d of Proposition 12.31 to the matrix M{ai, ...,«;). 

The second part follows also from Proposition 12.11 and from the obser- 
vation that the number of integral matrices satisfying the inequalities of 
Proposition 12.31 with d = N is finite. □ 

Corollary 2.5. A matrix of the form T{ai,a2, ■ ■ ■ ,cti) with Oi > 2 for 
i = 1,2, . . . ,1 is positive definite. 

Remark 2.6. Corollary \2. 51 follows of course also from the well-known ob- 
servation that T{2, 2,. . . ,2) is a Cram matrix for a root lattice of type A 
and from the fact that real symmetric positive matrices form a convex cone. 

Proof of Proposition 12. 3I Bv induction on I. The results hold obviously 
for / = 1 since 



We have obviously < b,d and d > b hy induction. Since < a < 6 we 
have 




satisfies all inequalites for ai >2. 
For the induction step we consider 




d 



b a — xb 
I —c + xd 



) 



a + xb> {x -l)b>b> {) 



for X > 2. 
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The inequalities < c < d show similarly 

-c + xd > (x - l)d > d > . 
The inequalities c — a < d — b and d > b imply 

{xd — c) — {xb — a) = x{d — b) — {c — a) > {x — l){d 
for X >2. 

3 Proof of Theorem 11.11 

Theorem 3.1. For every integral unimodular matrix 




such that 

< min(a, b, c, d) and max(5, c) < d 

there exists a unique integer I > 1 and a unique finite sequence (oi, . . . , ai) G 
{2,3,4, ...y such that 

A = M{ai)M{a2) ■ ■ ■ M{ai) . 
Theorem 11.11 is now implied by the following result. 

Proposition 3.2. Integral matrices of the form T{ai, . . . ,«„) with deter- 
minant N and diagonal coefficients > 2 are in bijection with the subset of 
integers in {1, . . . , N} which are invertible modulo N. 

Proof of Proposition 13.21 Given a pair (a. A*") of natural integers with 
a G {1,...,A^ — 1} such that a is invertible modulo A^, we consider b G 
{1, . . . , iV - 1} such that ab = 1 (mod N) and we set k = ^ € N. The 

matrix ^ = ( N ) hyptheses of Theorem 13.11 and has 

thus a unique factorisation of the form 

A = M{ai)M{a2) ■ ■ ■ M{ai) 

for some integer / > 1 and ai,a2, - ■ ■ G {2, 3, ... }^ 

Proposition 12.11 shows thus that the matrix T(ai, . . . , ai) has determi- 
nant N. 

This construction yields an injective map from (Z/A^Z)* into integral 
matrices of the form T{ai, a2, ■ ■ ■) with determinant and diagonal coef- 
ficients > 2. 



-b)>d-b>0 

□ 
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The fact that this map is onto fohows by applying the unicity result of 
Theorem 13.11 to the matrix 

-|r(a2,...,an-i)| -\T{a2,...,an)\\ _ ^j, . x 
\rrf M \rrt M ) — M[ai) ■ ■ ■ M{ai) 

associated to a suitable tridiagonal matrix T(ai, • • • )• '-' 
Proof of Theorem 13.11 If a = then unimodularity of A implies h = 
c = 1 and A = M{d) for some integer d'>2>\ = h = c. Otherwise, the 
identity 

he = ad + 1 

implies be = 1 (mod d) . This shows that the integers c, d are coprime and 
determine a, b uniquely by considering the unique integer b £ {1, . . . , d — 1} 
such that be = 1 (mod d) and by setting a = ^^j^ > 1. Unimodularity 
of A and strict positivity of a imply that Defining t € {1, ... ,6 — 1} as 
the unique integer satisfying the congruence t = —d (mod b), the equality 
det{A) = —ad + bd= 1 implies at = 1 (mod b) and we can thus consider 
the natural integer s = > 0. The computation 

-a -b \ f b t \ _ f Q -1 
c d J \ —a —s J \ 1 et — sd 

and the inequality 

b[et- , ^ d 1 = bet - adt + d = t + d > b 

show thus 



-a -b \ _ f -1 \ ( -t 
e d J \lei — sdj\a b 

where et — sd > 2. Since a = < b, the integral unimodular matrix 

-s -t 
a b 

satisfies again the hypotheses of Theorem 13. 11 Induction on d > 6 establishes 
the existence of a sequence (ai, a2, ■ ■ ■ ,ai) in {2, 3, 4, ... such that A = 
M{ai)M{a2)---M{ai). 

Unicity of the sequence ai, . . . is obvious for a = by Proposition 12.11 
For a > it follows from the observation that the inequalities < min(s, t) 
and t <b determine the integral unimodular matrix 



—s Oil ^ \ ( —a —b 

a b ~ [ -1 [ c d 



T{a2,a3,... 



uniquely. This shows ai = ai and by induction cti = and I' = I for 
ai, . . . , d;/ G {2, 3, . . . }' such that M(ai, . . . , a/) = M(di, . . . , d;'). □ 
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Remark 3.3. Theorem \3.1\ is equivalent to the assertion that the matrices 
M(2), M(3), M(4), . . . generate a free submonoid o/SL2(Z). The subgroup 
generated by M(2), M(3), ... is however not free and coincides with SL2(Z). 
Indeed, the identity 

M{\)M{\ + n)-^M{\) = M{\ - fi) 

shows that the subgroup generated by two integral matrices M{k) and M{k + 
1) contains the generators 

= ( ! "o^ ) = ( ? "i^ ) 

o/SL2(Z). 



4 Algorithmic aspects 

In the sequel, x (mod y) denotes always the unique natural integer in {1, . . . , y} 
representing the equivalence class of x in 'Ljy'L for an integer x and a nat- 
ural integer y >2. Similarly, if x is invertible modulo y, then x"^ (mod y) 
denotes the unique integer in {1, . . . ,y — 1} such that x{x~'^ (mod y)) = 1 
(mod y). 

Given a natural integer N > 2 and an integer a€ {!,... ,A^ — 1} which 
is invertible modulo A^, we denote by 

W{a,N) = ai,...,ai 

the unique finite sequence of natural integers > 2 such that 

M(ai...a;)= ( ~N ) ^ SL2(Z) 

where y = a^^ (mod N) and x = ""^^^ , cf Theorem 13.11 

We denote by l{a,N) = I the length of the word W{a,N) and by 
W{a, N)i = ai the i— th element (for i = 1, ... ,1) of the sequence W{a, N). 
We have thus 



Wia, N) = W{a, N)i,W{a, N)2,..., W{a, N\^^n) ■ 

We have obviously 

W{1,N) = W{l,N)i = N, lil,N) = 1 . 

Ifa, 6 G {!,..., are invertible elements modulo such that ab = 1 
(mod A^), Proposition 12.21 shows that we have 

l{a,N) = l{b,N) 
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and 

Otherwise stated, the matrices T{W {a^N)) and T{W(h,N)) are related to 
each other by conjugation with the antidiagonal involution of size /(a, N) x 
l{a,N). 

We denote by Jx[ the integral part of a real number x defined as the 
unique integer satisfying 

X — 1 < \x\ < X . 
The following result allows computations: 
Theorem 4.1. For every integer N > 2 we have 

1^(1, A^) = N . 

For an integer a € {2,... ,A^ — 1} invertible modulo N (with N a natural 
integer > 3) we have 





N 






a 



or equivalently 

W{a,N) 



1 + 



N 

T 



,W{{-N)-'^ {mod b),b) 



where b 



,-1 



(mod N)) E {!,..., iV- 1}. 



Remark 4.2. The map {a,N) i — > {{—N (mod a)), a) involved in the first 
equality for W{a,N) for a > 2 coincides (up to the — sign) with Euclid's 
celebrated algorithm for computing the greatest common divisor of a and N. 

Proof of Theorem 14.11 The assertion concerning W{1,N) is obvious. 
Set k = "^-^ > 1 where b = a^^ (mod N) is the inverse of a modulo N 
in {2, . . . , — 1}. Computing 



-k -b 
a N 



1 + [N/a\ 1 
-1 



-k{l + [N/a)\ +b -k 
a{l + [N/a\) - N a 



the lower left coefficient of the last matrix is given by 
N-{N (mod a)) 



a 1 



N = a+N-{N {mod a))-N = {-N) (mod a) 



and is thus a strictly positive natural integer < a. Since 

k{a{l + [N/a\) - N) = -kN = 1 (mod a) , 
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positivity of k,a{l + [A^/aJ) — N, the obvious inequality k < a and uni- 
modularity of all involved matrices imply that the last matrix satisfies the 
conditions of Theorem 14.11 This implies the first equality for a > 2 by 
induction on A^. 

The second equality can be proven similary. It follows also from the first 
equality and from Proposition 12.21 □ 

For computing W{a,N) with a G {1, . . . ,N — 1} invertible modulo N 
and very close to a huge integer the following result is useful: 

Proposition 4.3. The sequence W{N — l,N) is the constant sequence 
2,2, ... ,2 of length N - 1. 

We have for a G {2, . . . , — 2} invertible modulo N the equality 

W{a, N) = W{{X + l)a - AA^, Aa - (A - 1)N), 2, 2, 2, . . . , 2 

where A = la/{N — a) J and where 2, 2, 2, . . . , 2 is the constant sequence of 
length A. 

Proof The assertion for the sequence W{N — 1, A^) follows from the case 
k = N — 1 of the identity 

/ -1 y _ / l-zfc -k \ 
\l 2 J ~ \ k k+1 J 

which is easily established. 

For a > 2 such that a < N/2 we have a/ {N — o) < 1 yielding A = and 
the formula is trivial. 

The case a = N/2 implies {a,N) = (1,2) corresponding to the case 
{N -1,N) for N = 2. 

The proof in the case o > A^/2 is by induction on A. If A = 1 we have to 
prove 

W{a,N) = W{2a- N,a),2 

where ^N < a < ^N. This boils down to the identity l^(a, A^) = W{{-N) 
(mod a),a), (1 + [N/a]) of Theorem El 

For |A^ < a < A^ - 1 we have as above W{a, N) = W{2a - A^, a), 2 by 
Theorem 14.11 The trivial identities 

2a - A^ _ a 
a -{2a- N) ~ N -a ~ 

and 

k(2a -N)-{k- l)a = {k + l)a - kN 

and induction on A end the proof. □ 
Theorem 14.11 and Proposition 14.31 give a fast algorithm for computing 

l{a,N) 

E (W{a,N),r 
1=1 
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for arbitrary e S N (or e G C). The case e = corresponds to the length 
l{a,N) of the sequence W{a,N) and the case e = 1 yields the trace of the 
tridiagonal matrix T associated to W{a,N). The details are as follows: 

INPUT: (a, N) with a € {1, . . . , N -1} invertible modulo N (for iV > 2) 
and an integer e. 

Set r = 0. 

LOOP: 

IF a = 1, OUTPUT: (r + iV^) and STOP ENDIF. 

IF a < {N/2), replace r by r + (1 + [N/a\Y and {a,N) by {{-N) 
(mod a), a) where {—N) (mod a) is the choosen in {1, . . . , a — 1}) ENDIF. 

IF a = iV - 1, OUTPUT: r + {N - 1)2'= and STOP ENDIF. 

IF a > N/2, set A = la/{N - a)\, replace r by r + A2^ and {a,N) by 
((A + l)a - XN, Aa - (A - 1)A^) ENDIF. 

END OF LOOP. 



5 A function with an additional symmetry 

Set a{a,N) = Zk=i\w{a,N)k-l) = tiiT{W{a,N)))-dim{T{W{a,N))) 
for a S {1, . . . , N — 1} invertible modulo N. 

Proposition 12.21 shows that the function a satisfies a{a,N) = a{{a^^ 
(mod N)),N). 

The following result shows that a satifies an additional symmetry. 
Theorem 5.1. We have 

a{a,N) =a{N -a,N) 
for every a G {1, . . . , — 1} which is invertible modulo N. 

ProofWehavef7(l, Af) = N-l anda{N-l, N) = (iV-l)(2-l) = iV-1. 
Consider now a such that 2 < a < Setting A = [(A^ — a)/a] = 
i™"*^ ~ Proposition 14.31 shows 

aiN -a,N) = a{N - (A + l)a, N - Xa) + X 

implying 

a{N -a,N) = a{a, N - Xa) + X 
by induction on A^. Theorem 14.11 yields thus 

/,r // / 1 \ N N — Xa — (N (mod a)) 

a(N -a,N) =a((-N) (mod a), a + ^ ^ - + X 

a 

= a{{-N) (mod a),a) + [N/a\ 
which establishes the result since we have also 

a{a, N) = a{{-N) (mod a), a) + [A^/aJ 
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by Theorem 14.11 

A few values of cr(a, N) are given by the following Table: 



□ 



a = 
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10 


11 


12 


13 


N = 


2 : 


1 
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2 


2 
























N = 
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3 




3 






















N = 
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3 


4 




















N = 


6 : 


5 








5 


















N = 


7 : 
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6 
















N = 


8 : 


7 




4 




4 




7 














N = 


9 : 


8 


5 




5 


5 




5 


8 












N = 


10 


: 9 




5 








5 




9 










N = 


11 


: 10 
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5 
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6 


5 


5 


6 
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N = 


12 


: 11 








5 




5 








11 






N = 


13 


: 12 


7 


6 


6 


5 


7 


7 


5 


6 


6 


7 


12 




N = 


14 


: 13 




6 




6 








6 




6 




13 


N = 


15 


: 14 


8 




6 






8 


8 






6 




8 


N = 


16 


: 15 




7 




7 




6 




6 




7 




7 


N = 


17 


: 16 
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7 
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7 
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7 


6 


7 



16 



14 
7 



15 
9 
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5.1 Continued fraction expansions 

Let [71,72,...] be the continued fraction expansion of x £ (0,1) defined 
recursively by x = ^^^^.^ — j where 71 = [1/xJ and [72, ... ] = 1/x — 71 (and 
= [0] by convention). 

Theorem 5.2. We have 



l + aia,N) = Y,li 



i>l 



where 



[71 ) 72 5 • • • ] is the continued fraction expansion of the reduced 
rational fraction a/N in (0, 1). 

Proof The result holds obviously for a = 1. For a > 1, we have 

71 = [N/a\ 

and 

N-{N-{N (moda))) N (mod a) 

72, ... = N a - [N a\ = = 

a a 

The equality 71 = ai and induction on end the proof. □ 

Remark 5.3. Theorem \5.1\ follows also from Theorem \5.S\ together with the 
easy continued fraction expansion 

1 - X = [l,ai - 1,02,03, • • •] 

for x = [qi, a2! «3, • ■ • ] in {0, 1/2). 
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6 Lattices 



Given d integers ai, . . . , ^ {2, 3, ... j*^, we denote by A(ai, . . . , a^) the 
Euclideean lattice of rank d with scalar product given by the Gram matrix 
r(ai, . . . , Ud) (which is positive definite by Corollary [231 We say that such 
a lattice is of type T. 

Proposition 6.1. A lattice A{ai, . . . ,ad) of type T has minimal norm 
min(ai, . . . ,0^) > 2. 

Proof We work with a basis with Gram matrix T(ai, . . . , ad)- Replac- 
ing ai, . . . , by 2 yields a root lattice A' of type Ad and the obvious map 
sending a vector A = (Ai, . . . , A^) of the lattice A(ai, . . . , ad) to the corre- 
sponding vector A' = (Ai, . . . , A^) with the same coordinates in A(2, 2, . . . , 2) 
is a linear map which is one-to-one and length-shrinking. 

This proves the Proposition if min(ai, . . . , a^) = 2. In the remaining 
case, apply the Proposition to the lattice A' = A(qi — k, a2 — . . . , — k) 
where k = min(ai, . . . , a^) — 2 and observe that the obvious linear map 
from A(ai, . . . , ad) onto A' is one-to-one and diminuishes norms of non-zero 
elements at least by k. □ 

Proposition 6.2. Blocks of adjacents 2's in the sequence ai, . . . ,ad are in 

bijection with irreducible root-sublattices in A(ai, . . . , ad)- 

The easy proof is left to the reader. 

Theorem 6.3. Two lattices A(ai, . . . , ad) and A(&i, . . . , bd) of type T are 
isomorphic if and only if either ai = 6j for i = 1, . . . ,d or a^ = bd+i-i for 
i = 1, . . . , d. 

Corollary 6.4. The number of isomorphism classes of Euclidean lattices 
(of arbitrary dimension) of type T, determinant N > 2 and containing no 
elements of norm 1 is given by {^{N) + Vf^)/2 where equals the number 
of elements of order < 2 in the multiplicative group (Z/A^Z)*. 

Idea for the proof of Theorem 16.31 We show by induction on the 
dimension d that metric properties of a lattice A(ai, . . . , ad) of type T de- 
termine the integer sequence ai,...,ad uniquely, up to reversion of the 
order. 

If min(ai, . . . , ad) = 2, the root system TZ of A(ai , . . . , ad) is determined 
by the set of blocks of adjacent 2's in the sequence ai, 02, . . . • The sublat- 
tice A orthogonal to the root lattice generated by TZ corresponds to all terms 
ai such that ai-i,ai,ai+i > 2 and is an orthogonal sum of lattices of type 
T. By induction on the rank, the lattice A determines (up to reversion of 
the order) disjoint subsequences ai.,aij+i,ai.^d containing no 2's and not 
adjacent to 2's in ai, . . . , ad- Terms adjacent to irreducible root lattices can 
essentially be recovered by considering minimal norms of sets of vectors not 
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contained and not orthogonal to such an irreducible root lattice. Orthogo- 
nality consideration allow then to glue all pieces together in an essentially 
unique way. 

If min(ai, . . . , ctrf) > 2 one can replace the root lattice by the lattice 
generated by the set of minimal vectors (corresponding to basis vectors ztcj 
such that ai = min(ai, . . . , ad)) and proceed as above. □ 

Remark 6.5. Lattices of type T have poor densities and are thus not inter- 
esting from the point of view of sphere-packings. 

7 Polytopes 

Let A be an integral symmetric matrix of size d x d. For every natural 
integer N, we denote by Ca{N) the set of integral diagonal matrices of size 
dx d such that A + D \s positive definite of determinant N for every matrix 
D G Ca{N). We denote by Conv(/:^(iV)) C M'^ the convex hull of ^^(iV)- 

Theorem 7.1. For all N , the set Comr{CA{N)) is a polytope whose set of 
integral elements is contained in 

^T=NC-A{k) 

and whose vertices are given by the the integral elements Ca{N). 

Proof The Brunn- Minkowski Theorem, see eg. Theorem 6.2 in [3] states 
that 

det(^+ ^ XdD)^''^> Y1 Az5det(A + L»)i/'^ 

DeCA{N) DeCA{N) 

if ^^DeCAiN) = 1 with Ad positive real numbers. This shows that we 
have det{A-\-D') > N for all D' e Couv{Ca{N)) and the inequality is strict 
if D' ^ Ca{N). Finiteness of the set Ca{N) follows from the fact that given 
two elements (ai,...,a(i) and (/3i, . . . , /J^^) of Ca{N)., there exists indices 
i, j S {1, . . . , d} such that ai < f3i and aj > (3j. □ 

The construction of the polytope Conv(£/i(A^)) has the following prop- 
erties, given without proofs. (They are straightforward.) 

The diagonal part of A is essentially irrelevant: If A = A-\- D where D 
is integral diagonal and A is integral symmetric, then jC^{N) = Ca{N) — D. 
We suppose henceforth that A has zeroes along the diagonal. 

Automorphisms of A (conjugations by signed permutation matrices com- 
muting with A) act on the polytopes Conv(£^(A^)) in the obvious way by 
permuting the coefficients according to the underlying ordinary permutation 
matrix. 

If A is a "direct sum" A = Ai + A2 oi disjoint diagonal blocks, then 
Conv(£^(l)) = Conv(£^,(l)) x Conv(£^,(l)) . 
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More generally, 

Ca{N) = \Ji\nCaS) t^A.iN/l) 

where the union is over all natural integers divising N . 

The following construction yields a map from the permutation group 
on the d indices into the set Ca{^)'- Given a permutation matrix a of size 
d X d, we consider the unique integral diagonal matrix D = D-,^ such that 
the submatrix formed by the first k rows and columns of a~^{A + Dfj)a 
is unimodular for every integer k in {1, . . . These special vertices are 
invariant under automorphisms of A (and their convex hull is a polytope 
contained in Conv(>C^(l)) which is invariant under automorphisms of A). 

If A is an adjacency matrix of a simple finite undirected graph which is 
connected, the above construction can be restricted to permutations such 
that the first k vertices form connected subgraphs for all k. The subset of 
these vertices is still invariant under automorphisms. There exists slight 
generalisations of this construction. 

Examples The easiest case is given if the matrix A is the zero-matrix of 
size d X d. For d > 2, the polytope Conv(£^(A'')) is a (d — 1)— dimensional 
simplex if and only if is a prime number. More interesting cases occur if 
N is highly composite. 

A probably interesting example is given by considering the symmetric 
matrix A of size d x d with coefficients Aij = 1 if |i — j| = and Aij = 
otherwise. By [2], the set Ca{1) contains then exactly = + 1) 

elements forming the vertex set of an integral polytope Pd of dimension d if 
d >3. The polytope P3 for instance is the convex hull of the five vertices 

VI = (1,2,2) 
^^2 = (1,3,1) 
^^3 = (2,1,3) 
1^4 = (2,2,1) 
V5 = (3,1,2) 

Its faces are defined by equalities of the five inequalities 

y + z<4 {vi,V2,V3} 

x + y + z>5 {vi,V2,V4} 

x + 2y + z>7 {vi,V3,Vi,v^} 

2x + 3y + 2z<l3 {v2,V3,v^} 

x + y<4 {v2,va,v^} 

(with vertex-sets indicated for each face) and the polytope P3 is thus a pyra- 
mid with summit V2 over the planar polygone formed by the four vertices 
vi,V'i,V4^,v^. The face {vi,V2,v/^} contains all vertices associated to permu- 
tation matrices of size 3x3. (This can easily be generalised to an arbitrary 
dimension d: Such special vertices are all of the form (1, 2, 2, 2, . . . , 2, 1) -|- Cj 
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where = (0, . . . , 0, 1, 0, . . . , 0) denotes the i—th. basis vector and form a 
unique simphcial face of the of polytope Conv(£yi(l)).) 

It would perhaps be interesting to understand the polytopcs Pd in gen- 
eral, in particular eventual connections with so-called Associahedrons (or 
Stasheff polytopes) whose vertices are also enumerated by Catalan numbers 
and whose dimensions are in general one less. 
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